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ABSTRACT 

More than half of the atoms in the Universe recombined via forbidden transitions, 
so that accurate treatment of the forbidden channels is important in order to follow 
the cosmological recombination process with the level of precision required by future 
microwave anisotropy experiments. We perform a multi-level calculation of the re- 
combination of hydrogen (H) and helium (He) with the addition of the 2 3 P 1 -l 1 Sr J 
spin-forbidden transition for neutral helium (Hel), plus the nS-lS and nD-lS two- 
photon transitions for H (up to n = 40) and among singlet states of Hel (n ^ 10 
and 1^7). The potential importance of such transitions was first proposed by 
iDubrovich fc Grachev! (|2005l) using an effective three-level atom model. Here, we relax 
the thermal equilibrium assumption among the higher excited states to investigate the 
effect of these extra forbidden transitions on the ionization fraction x c and the Cos- 
mic Microwave Background (CMB) angular power spectrum Ce- The spin-forbidden 
transition brings more than a percent change in x c . The two-photon transitions may 
also give non-negligible effects, but currently accurate rates exist only for n ^ 3. We 
find that changes in both x e and Ci would be at about the percent level with the 
approximate rates given by IDubrovich fe Grachev! ((2005) . However, the two-photon 
rates from 3S to IS and 3D to IS of H appear to have been overestimated, and our 
best numerical calculation puts the effect on x c and CV at below the percent level. 
Nevertheless, we do not claim that we have the definite answer, since several issues re- 
main open; sub-percent level computation of the Cgs requires improved calculations of 
atomic transition rates as well as increasingly complex multi-level atom calculations. 

Key words: cosmology: cosmic microwave background - cosmology: early universe 
- cosmology: theory - atomic processes. 



1 INTRODUCTION 

The release of the third year data from the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) has further improved 
the precision with which we can constrain the cosmological 
parameters from the shape of the Cosmic Microwave Back- 
ground (CMB) anisotropies Ce i|Spergel et all [2 006). The 



Planck satellite, scheduled for launch in 2008 l|Planckll2006l) , 
will provide even higher precision Ce values and data down 
to smaller angular scales (£ < 2500). Higher precision in 
the observations requires increased accurarcy from the the- 
oretical calculations, in order for the correct cosmological 
parameters to be extracted. It now seems crucial to obtain 
the C'es down to at least the 1 percent level over a wide range 
oft 
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CMBFAST (|Seliak fc Zaldarriagalll996l ) is the most com- 
monly used Boltzmann code for calculating the Ces, and 
it gives consistent results with other independent codes (see 
Seliak et al.ll200 and references therein) . The dominant un- 
certainty in obtaining accurate Ces comes from details in the 
physics of recombinat ion, for example, the 'fudge factor' i n 
the recfast routine l|Seager. Sasselov fc Scottll 19991 , |2000|) . 
Calculatio ns of co s molog ical r ecombination were first pub- 
lished by Peebles! (Il968 | ) and Zel'dovich. Kurt fc Sunvaevl 



|l968h . ISeager. Sasselov fc Scott (2000) presented the most 
detailed multi-level calculation and introduced a fudge fac- 
tor to reproduce the results within an effective three-level 
model. Although the multi-level calculation already gives 
reasonable accuracy, the required level of accuracy contin- 
ues to increase, so that today any effect which is ~ 1 per cent 
over a range of multipoles is potentially significant. Several 
modifications have been recently suggested to give per cent 
level changes in the ionization fraction and/or the Ces (see 
Section 4 for details). Most of these modifications have been 
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calculated only with an effective three-level code, and so 
the results may be different in the multi-level calculation, 
since there is no thermal equilibrium assumed between the 
upper states. Here we want to focus on one of these mod- 
ific ations, namely the ext r a forb idden transitions proposed 
bv iDubrovich fc Grachev! ([Ill), which we study using a 
multi-level code. 

In the standard calculations of recombination, one 
considers all the resonant transitions, but only one for- 
bidden transition, which is the 2S-1S two-photon tran- 
sition, and this can be in cluded for both H and He. 
IDubrovich fc Grachev! (|2005h suggested that one should also 
include the two-photon transitions from higher excited S 
and D states to the ground state for H and He I, and also 
the spin-forbidden transition between the triplet 2 3 Pi and 
singlet ground state 1 So for He I. They showed that the 
recombination of both H and He I sped up in the three- 
level atom model. The suggested level of change is large 
enough t o bias the determination of t he cosmological pa- 
rameters l|Lewis. Weller fc Battvell2006h . 

The purpose of this paper is to investigate the 
effe ct of the extra forbidde n transitions suggested 
by IDubrovich fc Grachev! (|2005h in the multi-level atom 
model without assuming thermal equilibrium among the 
higher excited states. The outline of the paper is as follows. 
In Section 2 we will describe details of the rate equations in 
our numerical model. In Section 3 we will present results on 
the ionization fraction x c and the anisotropies Ce, and assess 
the importance of the addition of the forbidden transitions. 
Other possible improvements of the recombination calcula- 
tion will be discussed in Section 4. And finally in Section 5 
we will present our conclusions. 



2 MODEL 

In this paper we follo w the formalism of the mu l ti-leve l cal- 
culation performed bv lSeager. Sasselov fc Scottl l|2000l ). We 
consider 100 levels for HI, 103 levels for He I, 10 levels for 
Hell, 1 level for He III, 1 level for the electrons and 1 level 
for the protons. For H, we only consider discrete n levels 
and assume that the angular sub-levels (^-states) are in sta- 
tistical equilibrium within a given shell. For both He I and 
Hell, we consider all the i'-states separately. The multilevel 
He I atom includes all states with n 10 and £ ^ 7. Here we 
give a summary of the rate equations for the number den- 
sity of each energy level i, and the equation for the change 
of matter temperature Tm- The rate equation for each state 
with respect to redshift z is 

(1 + *)^= 1 
az 



H{z) 

N 

(n c n c R C i — ruRic) + Yj ARji 



+ 3n,, (1) 



where rn is the number density of the ith excited atomic 
state, n e is the number density of electrons, and n c is 
the number density of continuum particles such as a pro- 
ton, Hell, or He III ion. Additionally R C i is the photo- 
recombination rate, Ri c is the photo-ionization rate, ARji 
is the net bound-bound rate for each line transition, and 
H(z) is the Hubble parameter. We do not include the 



collisi onal rates, as they have be en shown to be negligi- 
ble (|Seager, Sasselov fc Scotll2000h . 

Fo r He I, we update the atom ic data for the energy 
levels (jMorton. Wu fc Drake! l2006h. the os cillator strength 
for resonant transitions ( Drake fc Mortonl . in preparation) 
and the photo-ionization cross-section sp ectrum. We use the 
hoto-ionization cross-section given by iHummer fc Storevl 
1998) for n ^ 10 and ^^4, and adopt the hydrogenic ap - 
proximation for states with £ ^ 5 (jStorev fc Hummerlll99ll ) . 
It is hard to find published accurate and complete data 
for the photo-ionization cross-section of He I with large n 
and I. For example, a recent paper bv lBauman et all (|2005t) 
claimed that they had calculated the photo-ionization cross- 
section up to n = 27 and £ = 26, although, no numerical val- 
ues were provided. 

The rate of change of matter temperature with respect 
to redshift is 



(1+2! 



dz 



8a T U 



3H(z)m e c n c + n H + n H , 



■(Tm-T r )+2T m ,(2) 



where Tr is the radiation temperature, riHe is the total num- 
ber density for helium, m e is the electron mass, c is the speed 
of light, U — arTr, or is the radiation constant and ctt is 
the Thompson scattering cross-secti on. 

ISeager. Sasselov fc Scottl (|2000h considered all the res- 
onant transitions and only one forbidden transition, namely 
the 2S-1S two-photon transition, in the calculation of each 
atom, (for Hel, 2S = 2 1 S and IS = l x So). The 2S-1S two- 
photon transition rate is given by 



AR 2 



A 2S ( n 2S -nis-e-'' p " 2S - ls/tBTM 
Sis 



(3) 



where A2S is the spontaneous rate of the corresponding two- 
photon transition, i^s-is is the frequency between levels 2S 
and IS, gi is the degeneracy of the energy level i, hp is 
Planck's constant and &b is Boltzmann's constant. 

Here we include the following extra forbidden transi- 
tions, which were first suggested by IDubrovich fc Grachev! 
(2005). The first ones are the two-photon transitions from 
nS and nD to IS for H, plus n 1 So and n 1 T>2 to l 1 So for Hel. 
For example, for H, we can group together the nS and nD 
states coming from the same level, so that we can write the 
two-photon transition rate as 

A-RnS+nD^lS = AjJs + nD X 



SnS + nD -h r v nl /k B T N 
HnS+nD - niS e 

Sis 



■(4) 



Here n (without a subscript) is the principle quantum num- 
ber of the state, n n s+„D is the total number density of the 
excited atoms in either the nS or nD states, and A nS+ „ D is 
the effective spontaneous rate of the two-photon transition 
from nS + nD to IS, which is approxi mated by the following 
formula (|Dubrovich fc Grachev! I2005T ): 



A 



nS + nD 



54A? S 



SnS + nD V n + 1 



lln 1 - 41 



(5) 



where A^ R is equal to 8.229 s" 1 (|Goldmanl Il989l : 
ISantos. Parente fc Indelicate! 19981). The latest value of AS 



is equal to 8.2206 s 1 1 Labzowskv. Shonin fc Solovve v booa i 
and does not bring any noticeable change to the result. 
Here <?„s+nD is equal to 1 for n = 2, and 6 for n ^ 3. This 
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spontaneous rate is estimated by considering only the non- 
resonant two-photon transition s thro ugh one intermediate 
state nP. iDubrovich fc Grachev! (|2005h ignored the resonant 
two-photon transition contributions, since the escape prob- 
ability of these emitted photons is very low. The above for- 
mula for AjJ s+nD is valid up to n ~ 40, due to the dipole 
approximation used, although it is not trivial to check how 
good this approximate rate is. Besides the 2S-1S two-photon 
rate, only the non-resonant two-photon rates from 3S to IS 
and 3D to IS are calcula t ed ac curately and available in the 
literature. ICresser et al.l (|l986l ) evaluated Ajj? s and A^ D by 
including the non-resonant transitions through the higher- 
lying intermediate nP states (ra^4), which are equal to 
8.2197 s _1 an d 0.13171s" 1 , respectively. These values were 
confirmed by iFlorescu. Schneider fc Mihailescul (|l988l ) and 
agreed to three significant figures. Using these values, we 
find that A^ s+nD is equal to 1.484 s -1 , which is an order 
of magnitude smaller than the value from the approximated 
rat e coming from equation (5). The approximation given 
by IDubrovich fc Grachevl {20051 ) therefore seems to be an 
overestimate. This leads us instead to consider a scaled rate 
AjJs+nD i which is equal to A^ s+nD multiplied by a factor to 
bring the approximated two-photon rates of H (equation (5)) 
with n = 3 into agreement with the numerical value given 
above, i.e. 



A„ s+nD = 0.0664 A 



(6) 



Note that the use of the non-resonant rates is an ap- 
proximation. The resonant contributions are suppressed in 
practice because of optical depth effects, and in a sense some 
of these contributions are already included in our multi-level 
calculation. Nevertheless, the correct way to treat these ef- 
fects would be in a full radiative transfer calculation, which 
we leave for a future study. For He I, we treat n 1 So and 
n D 2 separately and use equation (3) for calculating the 
tr ansition rates. The spon taneous rate A„|y nD is estimated 
by IDubrovich fc Grachevl 1120051 ) by assuming a similar form 
to that used for A^ s+nD : 



A 



S/nD 



1045yl Hcl 



n+1 



lln 2 



41 



(7) 



where A HcI is a fitting parameter (which is still uncer- 
tain both theoretically and experimentally). According to 
IDubrovich fc Grachevl {2005), resonable values of A range 
from 10 to 12s _1 , and we take A = lis -1 here. In our calcu- 
lation, we include these extra two-photon rates up to n — 40 
for H and up to n = 10 for He I. 

The other additional channel included is the spin- 
forbidden transition between the triplet 2 3 Pi and sin- 
glet l 1 So states in Hel. This is an intercombination/semi- 
forbidden electric-dipole transition which emits a single pho- 
ton and therefore we can calculate the corresponding net 
rate by using the bound-bound resonant rate expression, 
i.e. 

A_R 2 3 Pl _ii So = p 8 3p 1 ,iis x 

(n 2 3 Pl ii 2 3p lil i So - n 2 i So _R 1 i Soj2 3 Pl ) , (8) 



where 



^Pi.^So = yl 2 3 P 1 ,i 1 s + ^Pi.^So^ 



(9) 
(10) 



Figure 1. The ionization fraction x c as a function of redshift z. 
The solid line is calculat e d usin g the original multi-level code of 
ISeager. Sass clov & Scott (2000), while the dashed line includes 



all the extra forbidden transitions discussed here. 



^2 3 Pi,l 1 S ^23Pi,l 1 S 
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ff23p t 



Si 1 S 



"23 P X 



(11) 

(12) 



Here A 9 



2 3 p 1 ,i 1 s > - B 2 3 p 1 ,i 1 s and £ii Soi2 3 Pl 
coefficients, p 2 3 Pll i So is the Sobo lev escape probability, r a 
is th e Sobolev optical depth (see ISeager. Sasselov fc Scott] 
20001 . and references therein), A 2 3p lil ig is the wavelength 
of the energy difference between states 2 3 Pi and l 1 So, and 
J is the blackbody intensity with temperature Tr. 

This 2 3 Pi-l 1 So transition is not the lowest tran- 
sition between the singlet and the triplet states. The 
lowest one is the forbidden magnetic-dipole transi- 
between 2 3 Si and l 1 So, with Einstein coefficient 
1.73xl0 _4 s _1 {Lin. Johnson fc Dalgarnol 
However, this 



are the Einstein 



tion 



■A a 3H .,]lS 



.4 



2 3 P 1 .1 l Sn 



177.58 s" 



i s much sma ller than 
{Lach fc Pachuckil 1200 ll ; 
preparation) , so this transition 
Note that IDubrovich fc Grachevl 



Drake fc Morton! m 
can be neglected. 

20051) used an older value of j4 2 3 Pi 1 i So = 233 s 



Lin. Johnson fc Dalgarnolll977l ) in their calculation. 



We use the Bader-Deuflhard semi- implicit numerical in- 
tegration scheme (see Section 16.6 in iPress et al.l fl992h to 
solve the above rate equations. All the numerical results 
are carried out using the ACDM model with cosmologi- 
cal parameters: f2p=0.04; J7cdm = 0.2; S7a = 0.76; S7k = 0; 
y p =0.24; Tp = 2-725 K and h = 0.73 (consistent with those 
m ISpergel et all kooe). Here Y p is the primordial He abun- 
dance and To the present background temperature. 



3 RESULTS 

3.1 Change in ionization fraction 

The recombinatio n histories calculated using the previous 
multi-level code {Seager. Sasselov fc Scott] |2000| ) and the 
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Figure 2. The fractional difference ('new' minus 'old') in x e be- 
tween the two models plotted in Fig. [T] as a function of redshift 
2. The solid and dot ted lines are the model s with the two- photon 
rates for H given bv lDubrovich fc Grachevl l|2005t) and the scaled 
one given by equation (6), respectively. Both curves are calculated 
using all the He I forbidden transitions as discussed in the text. 

code in this paper are shown in Fig. [T] where x c = n e /nn 
is the ionization fraction relative to hydrogen. As we have 
included more transitions in our model, and these give elec- 
trons more channels to cascade down to the ground state, 
we expect the overall recombination rate to speed up, and 
that this will be noticeable if the rates of the extra frobidden 
transitions are significant. From Fig. [1] we can see that the 
recombination to He I is discernibly faster in the new calcu- 
lation. Fig. [5]shows the difference in x c with and without the 
extra forbidden transitions. The dip at around z = 1800 cor- 
responds to the recombination of He I and the one around 
z — 1200 is for H. Ov erall, the addition of the for bidden 
transitions claimed by iDubrovich fc Grachevl (1200 5l ) leads 
to greater than 1 per cent change in x c over the redshift 
range where the CMB photons are last scattering. 

In the last Section, we found that the approxim ated 
two-photon rate given bv lDubrovich fc Grachevl (|2005l ) for 
H with n = 3 was overestimated by more than a factor of 10. 
By considering only this extra two-photon transition, the 
approximate rate gives more than a per cent difference in 
x e , while with the more accurate numerical rates, the change 
in x c is less than 0.1 per cent (as shown in Fig. |3}. Based 
on this result, we do not need to include this two-photon 
transition, as it brings much less than a per cent effect on x e . 
For estimating the effect of the extra two-photon transitions 
for higher n, we use the scaled two-photon rate given by 
equation (6). The result is plotted in Fig. [4] The change in 
a; e with the scaled two-photo n rates is no more than .4 per 
cent, while the one with the IDubrovich fc Grachevl (|2005l ) 
approximated rates brings about a 5 p er cen t change. 

For He I, IDubrovich fc Grachevl (|2005h included the 
two-photon transitions from n = 6 to 40, since they claimed 
that the approximate formula (equation 7) is good for n > 6. 
In our calculation, we use A^ nD from the approximate for- 
mula for the two-photon transitions of n = 3 to 10, since this 
is the best one can do for now (and the formula at least gives 
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Figure 3. Fractional change in x e with the addition of the two- 
photon transition from 3S and 3D to IS for H. The solid line is cal- 
culate d with the approximate rate given bvlDubrovich ~fc Grachevl 
l|2005h while the dashed line is calculated with the numerical rates 
given bv lCresser et al" 

the right order of magnitude). The addition of the singlet- 
triplet 2 3 Pi-l 1 So transition and the n 1 So-l 1 So and n 1 T>2~ 
l 1 So two-photon transitions with n = 3—10 cause more than 
1 per cent changes in x c (as shown in Fig. 2). The 2 3 Pi-l 1 So 
transition has the biggest effect on x D . 

Fig. [5] shows the fractional difference in x e using dif- 
ferent combinations of additional forbidden transitions. We 
can see that the 2 3 Pi-l 1 So transition alone causes more 
than a 1 per cent change in x c , and the addition of each 
two-photon transition only gives about another 0.1 per cent 
change. The extra two-photon transitions from higher ex- 
cited states (larger n) have a lower effect on x e compared 
with that from small n, and we checked that this trend con- 
tinues to higher n. However, the convergence is slow with 
increasing n. Therefore, one should also consider these two- 
photon transitions with n > 10 for He I, and the precise result 
will require the use of accurate rates, rather than an ap- 
proximate formula such as equation (7). For the 2 3 Pi-l 1 So 
transition, IDubrovich fc Grachevl (2005) adopted an older 
and slightly larger rate, and this causes a larger change of 
the ionization fraction (about 0.5 per cent more compared 
with that calculated with our best rate), as shown in Fig. [5] 



3.2 The importance of the forbidden transitions 

One might wonder why the semi-forbidden transitions are 
significant in recombination at all, since the spontaneous 
rate (or the Einstein A coefficient) of the semi-forbidden 
transitions are about 6 orders of magnitude (a factor of a 2 , 
where a is the fine-sturcture constant) smaller than those 
of the resonant transitions. Let us take He I as an example 
for explaining the importance of the spin-forbidden 2 3 Pi- 
l 1 So transition in recombination. The spontaneous rate is 
equal to 177.58 s~ for this semi- forbidden transition, which 
is much smaller than 1.7989 x 10 9 s _1 for the 2 1 Pi-l 1 S 
resonant transition. But when we calculate the net rate 
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Table 1. The percentage of electrons cascading down in each channel from n = 2 states to 
the 1 So ground state for He I. 



2 1 S -> 1% 


2!Pi -► l^o 


2 3 Pi -> liSo 


(two-photon) 


(resonant) 


(spin-forbidden) 


code from Seaeer, Sasselov & Scott (2000) 30.9% 


69.1% 




this work 17.3% 


39.9% 


42.8% 



-ft 




Our best value 

/ Dubrovich & Grachev 



Figure 4. Fractional change in x c with the addition of different 
forbidden transitions for H. The long-dashed, dotted, dashed and 
solid lines include the two-photon transitions up to n = 10, 20, 30 
and 40, respectively, using the approximation for the rates given 
by equation (5). The dot-dashed line is calculated with the scaled 
rate from equation (6). 
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Figure 5. Fractional change in x e with the addition of different 
forbidden transitions for Hel as a function of redshift. The solid 
line corresponds to the calculation with only the 2 3 Pi— l^o spin- 
forbidden transition. The short-dashed, dotted, long-dashed, dot- 
dashed and long dot-dashed lines include both the spin-forbidden 
transition and the two-photon (27) transition(s) up to n = 3, 4, 
6, 8 and 10, respectively. 



Figure 6. Fractional change in x e with only the 2 3 Pi— ^So 
forbidden transition. T he solid line is calcu lated with the 
rate A = 233s — 1 from iDubrovich fc Grachev! |200d ) and the 
dash ed line is computed wi th our best value A = 177.58s — 1 
from lLach fc Pachuck] j200lF ) . 




Figure 7. Escape probability pij as a function of redshift. The 
solid line corresponds to the resonant transition between 2 1 ?! and 
l^o, while the dashed line refers to the spin-forbidden transition 
between 2 3 Pi and l 1 Sn. 
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Figure 8. Net bound-bound rates for He I as a function of red- 
shift. The solid line is the resonant transition between 2*Pi and 
l^So, the short-dashed line is the two-photon transition between 
2 1 Si and l^So- And the long-dashed line is the spin-forbidden 
transition between 2 3 Pi and 1 1 Sq. 



[see equation (8)], we also need to include the effect of ab- 
sorption of the emitted photons by the surrounding neu- 
tral atoms, and we take this into account by multiplying 
the net bound-bound rate by the So bolev escape probability 
Pij (|Seager. Sasselov fc Scottl 2000). If pij — 1, the emitted 
line photons can escape to infinity, while if Pij = the pho- 
tons will all be reabsorbed and the line is optically thick. 
Fig. [7] shows that the escape probability of the 2 1 Pi-l 1 So 
resonant transition is about 7 orders of magnitude smaller 
than the spin-forbidden transition. This makes the two net 
rates roughly comparable, as shown in Fig. [5] From equa- 
tion (11), we can see that the easier it is to emit a photon, the 
easier that photon can be re-absorbed, because the optical 
depth t s is directly proportional to the Einstein A coefficient. 
So when radiative effects dominate, it is actually natural to 
expect that some forbidden transitions might be important 
(although this is not true in a regime where collisonal rates 
dominate which is often the case in astrophysics) . In fact for 
today's standard cosmological model, slightly more than half 
of all the hydrogen a toms in the Universe recomb ined via a 
forbidden transition l|Wong. Seager fc Scottl l200fj ). Table 1 
shows that this is also true for helium. 

Iri the previous multi-level calcula- 
tion (Scagcr, Sass elov fc Scottl l2000h . there was no direct 
transition between the singlet and triplet states. The only 
communication between them was via the continuum, 
through the photo-ionization and photo-recombination 
transitions. Table 1 shows how many electrons cascade 
down through each channel from n = 2 states to the ground 
state. In the previous calculation, about 70% of the electrons 
went down through the 2 1 Pi-l 1 So resonant transition. In 
the new calculation, including the spin-forbidden transition 
between the triplets and singlets, there are approximately 
the same fraction of electrons going from the 2 1 Pi and 2 3 Pi 
states to the ground state (actually slightly more going 
from 2 3 Pi in the current cosmological model). This shows 
that we should certainly include this forbidden transition in 
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Figure 9. Relative change in the temperature (TT) angular 
power spectrum due to the addition of the forbidden transitions. 
The solid line includes the spin-forbidden transition and also the 
two-photons transitions up to n = 10 for He I, the dotted line 
includes all the above transitions and also the two-photon tran- 
sitions u p to n = 40 for H calcu l ated w ith the approximate rates 
given by iDubrovich fc Grachevl d2005r i . The dashed line is com- 
puted with the same forbidden transitions as the dotted line, but 
with our scaled rates (and represents our best current estimate). 



future calculations. Our estimate is that only about 40% of 
helium atoms reach the ground state without going through 
a forbidden transition. 

How about the effect of other forbidden transitions 
in He I recombination? We have included all the semi- 
forbidden electric-dipole transitions with n ^ 10 and £ ^ 7, 
and with oscilla tor strengths larger than 10 -6 given by 
iDrake fc Morton] (private communication). There is no sig- 
nificant change found in the ionization fraction. Besides 
the 2 3 Pi-l 1 So transition, all the other extra semi-forbidden 
transitions are among the higher excited states where the 
resonant transitions dominate. This is because these transi- 
tion lines are optically thin and the escape probabilities are 
close to 1. 



3.3 Effects on the anisotropy power spectrum 

The CMB anisotropy power spectrum Ce depends on 
the detailed profile of the evolution of the ionization 
fraction x c . This determines the thickness of the pho- 
ton last scattering surface, through the visibility function 
g(z) = e~ T dr/dz, where r is the Thomson scattering opti- 
cal depth (r = cot / n c (dt / dz) dz) . The function x a (z) sets 
the epoch when the tight coupling between baryons and pho- 
tons breaks down, i.e. when the photon diffusion length be- 
comes long, and the visibility function fixes the time when 
the fluctuations are effectively frozen in (see lHu et alJll995l : 
ISeager. Sasselov fc S cott 2000, and references therein). The 
addition of the extra forbidden transitions speeds up both 
the recombination of H and He I, and hence we expect that 
there will be changes in Ci. 

In order to perform t he required calculation, we have 
used the code cmbfast l|Seliak fc Zaldarriagal Il996t ) and 
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Figure 10. Relative change in the polarization (-E-E) angular 
power spectrum due to the addition of the forbidden transitions, 
with the curves the same as in Fig. [9] 



modified it to allow the input of an arbitrary recombina- 
tion history. Figs. l9l and 1 101 show the relative changes in the 
CMB temperature (TT) and polarization (EE) anisotropy 
spectra, respectively, with different combinations of extra 
forbidden transitions. The overall decrease of free electrons 
brings a suppression of Ce over a wide range of I. 

For He I, there is less x c at z ~ 1400 — 2500, which leads 
to an ealier relaxation of tight coupling. Therefore, both the 
photon mean free path and the diffusion length are longer. 
Moreover, the decrease of x c in the high-z tail results in in- 
creased damping, since the effective damping scale is an av- 
erage over the visibility function. This larger damping scale 
leads to suppression of the high-i? part of the power spec- 
trum. From Figs. 151 and mil we can see a decrease of Ce (for 
both TT and EE) toward high I for He I, with the maximum 
change being about 0.6 percent. 

For H, the change of Ct is due to the decrease in x c at 
z~600 — 1400 (see Fig. [2}. There are two basic features in 
the curve of change in Ce (the dotted and dashed lines in 
Fig. |5J). Firstly, the power spectrum is suppressed with in- 
creasing £, due to the lower x e in the high- 2 tail (z > 1000). 
Secondly, there are a series of wiggles, showing that the lo- 
cations of the acoustic peaks are slightly shifted. This is due 
to the change in the time of generation of the Ces in the 
low- z tail. Cf E actually shows an increase for t ^ 1000 (see 
Fig. HOP : this is caused by the shift of the center of the visibil- 
ity function to higher z, leading to a longer diffusion length. 
Polarization occurs when the anisotropic hot and cold pho- 
tons are scattered by the electrons. The hot and cold photons 
can interact with each other within the diffusion length, and 
therefore, a longer diffusion length allows more scatterings 
and leads to a higher intensity of polarization at large scales. 

With the app roxim ate rates used 

by iDubrovich fc Grachevf {2005), the maximum rela- 
tive change of Cf T is about 4 percent and for Cf E it is 
about 6 percent. The overall change is thus more than 1 
percent over a wide range of I. However, if we adopt the 
scaled two- photon rate given by equation (6), the relative 
changes of Cj T and Cf are no more than 1 per cent. Note 



that we do not plot the temperature-polarization correla- 
tion power spectrum here, since there is no dramatically 
different change found (and relative differences are less 
meaningful since Cj E oscillates around zero). 



4 DISCUSSIONS 

In our model we only consider the semi-forbidden transitions 
with n ^ 10 and I ^ 7 for He I and the two-photon transitions 
from the higher S and D states to the ground state for H 
and He I. It would be desirable to perform a more detailed 
investigation of all the other forbidden transitions, which 
may provide more paths for the electrons to cascade down 
to the ground state and speed up the recombination process. 
In this paper we have tried to focus on the forbidden tran- 
sitions which are likely to be the most significant. However 
we caution that, if the approximations used are inadequate, 
or other transitions prove to be important, then our results 
will not be accurate. 

There are several other approximations that we have 
adopted in order to perform our calculations. For example, 
we consider the non-resonant two-photon rates for higher 
excited rates. The two-photon transitions from higher ex- 
cited states (n ^ 3) to the ground state are more compli- 
cated than the 2S-1S transition, because of the resonant 
intermediate states. For example, for the 3S-1S two-photon 
transition, the spectral distribution of the emitted photons 
shows infinities (resonance peaks) at the freque ncies corre- 
spon ding to the 3S-2P and 2P-1S transitions (|Tung et al.l 
1 19841 ). Here, we use only the non-resonant rates, by assum- 
ing a smooth spectral distribution of the emitted photons; 
this probably gives a lower limit on the change of x c and Ce 
coming from these extra forbidden transitions. The correct 
way to treat this would be to consider the rates and feed- 
back from medium using the full spectral distribution of the 
photons and radiative transfer; this will have to wait for a 
future study. 

Besides the consideration of more forbidden tran- 
sitions, there are many other improvements that could 
be made to the recombination cal c ulatio n. In particular, 
iRubino-Martin. Chluba fc Sunvaevl ((2006) showed that a 
multi-level calculation of the recombination of H with the 
inclusion of seperate ^-states can give more than 20 per cent 
difference in the population of some levels compared with 
the thermal equilibrium assumption for each n-shell. The 
lat est calculation, considering up to 100 shel ls, is presented 
bv lChluba. Rubino-Martin fe Sumrae^l (|2006l ). but does not 
include all the forbidden transitions studied here. A more 
complete calculation should be done by combining the for- 
bidden transitions in a code with full angular momentum 
states, and we leave this to a future study. There are also 
other elaborations which could be included in future calcu- 
lations, which we now describe briefly. 

The rate equation we use for all the two-photon tran- 
sitions only includes the spontaneous term, assuming there 
is no inter action with the radi a tion b ackground (see equa- 
tion (3)). I Chluba fc Sunvaevl (|2006l ) suggested that one 
should also consider the stimulated effect of the 2S-1S two- 
photon transition for H, due to photon s in the low frequency 
tail of the CMB blackbody spectrum. iLeung. Chan fc Chul 
l|2004l ) additionally argued that the change of the adiabatic 
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index of the matter should also be included, arising due to 
the neutralization of the ionized gas. These two modifica- 
tions have been studied only in an effective three-level atom 
model, and more tha n a percent cha n ge in x e was claimed 
in each case (but see IWong fc Scott) |2006t ) for arguments 
against the effect claimed by Leung et al. 2004). 

For the background radiation field J, we approximated 
it with a perfect blackbody Planck spectrum. This approx- 
imation is not completely correct for the recombination of 
H, since the He line distortion photons redshift into a fre- 
qu ency range that can in principle photo-ionize the neutral 
H dDell' Antonio fc Rvbickil Il 993: Scag er. Sasselov fc Scott) 
l2000l : I Wong Seaeer fc ScottJ 120061 ) . Althought we expect 
this secondary distortion effect to bring the smallest change 
on x a among all the modifications suggested here, it is 
nevertheless important to carry out the calculation self- 
consistently, particularly for the spectral line distortions. 
In order to obtain an accurate recombination history, we 
therefore need to perform a full multi-level calculation with 
seperate ^-states and all the improvements suggested above, 
which we plan to do in a later paper. 

For completeness we also point out that the accu- 
racy of the physical constants is important for recom- 
bination as well. The most uncertain physical quantity 
in the recombination calculation is the gravitational con- 
stant G. The value of G used previously in the REC- 
fast code is 6.67259 x 1 0~ 11 m 3 kg~ 1 s~ 2 and the la t est va lue 
(e.g. from the Particle Data Group, lYao et all l2006h is 
6.6742 x I0 _11 m 3 kg~ 1 s -2 . Another quantity we need to 
modify is the atomic mass ratio of 4 He and 1 H, m4 He /mi H , 
wh ich was previou sly taken to be equal to 4 (as pointed out 
bvlSteigmanll20 06). By using the atomic masses given by 
lYao et aL 1 20061 ). the mass ratio is equal to 3.9715. The 
overall change in x c is no more than 0.1 per cent after up- 
dating these two constants in both recfast and multi-level 
code. 



5 CONCLUSIONS 

In this paper, we have computed the cosmological recombi- 
nation history by using a multi-level code with the addition 
of the 2 3 Pi to l 1 So spin-forbidden transition for He I and 
the two-photon transitions from nS and nD states to the 
ground st ate for both H and He I. Wi th the approximate 
rates from iDubrovich fc Grachev! |2005), we find that there 
is more than a per cent decrease in the ionization fraction, 
which agrees broadly with the result they claimed. How- 
ever, the only available accurate numerical value of two- 
photon rate with n ^ 3 is for the 3S to IS and 3D to 
IS t ransitons for H. We found th at the approximate rates 
from IDubrovich fc Grachev! (|2005l ) were overestimated, and 
instead we considered a scaled rate in order to agree with 
the numerical n = 3 two-photon rate. With this scaled rate, 
the change in x c is no more than 0.5 per cent. 

Including these extra forbidden transitions, the change 
in the CMB anisotropy power spectrum is more than 1 per 
cent, which will potentially affect the determination of cos- 
mological parameters in future CMB experiments. Since one 
would like the level of theoretical uncertainty to be negli- 
gible, it is essential to include these forbidden transitions 
in the recombination calculation. In addition, we still re- 



quire accurate spontaneous rates to be calculated for the 
two-photon transitions and also a code which includes at 
least all the modifications suggested in Section 4, in order 
to obtain the Cts down to the 1 per cent level. Achieving 
sub-percent accuracy in the calculations is challenging! 

However, the stakes are high - the determination of the 
parameters which describe the entire Universe - and so fur- 
ther work will be necessary. Systematic deviations of the 
sort we have shown would potentially lead to incorrect val- 
ues for the spectral tilt derived from Planck and even more 
ambitions future CMB experiments, and hence incorrect in- 
ferences about the physics which produced the density per- 
turbations in the very early Universe. It is amusing that in 
order to understand physics at the 10 15 GeV energy scale we 
need to understand eV scale physics in exquisite detail! 
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